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ABSTRACT. This paper reviews selected aspects of the growth of cosmologi- 
cal structure, covering the following general areas: (1) expected characteristics 
of linear density perturbations according to various candidate theories for the 
origin of structure; (2) low-order theory for statistical measures of fluctua- 
tions; (3) formation of nonlinear structures and nonlinear evolution of the 
mass distribution; (4) the relation between the density field and the galaxy 
distribution; (5) constraints on cosmological models from galaxy clustering 
and its evolution. 



1 Background 

This conference takes place at a 
time when the subject of large-scale 
structure is reaching a certain matu- 
rity, following a period of very rapid 
development. To appreciate how far 
we have come, it is instructive to look 
back nearly 20 years to the 1979 Les 
Houches summer school on Physical 
Cosmology, whose proceedings were 
greatly influential in their time. 

In 1979, there was an apprecia- 
tion that galaxies displayed roughly 
power-law correlations (measured on 
small scales only), and that these 
would only grow by gravity if there 
were some significant form of initial 
fluctuations. However, the central 
problem of the origin of these initial 
fluctuations was scarcely mentioned. 
The only candidate primordial spec- 
tra mentioned were expressed in 
terms of the density fluctuation in a 
region containing n particles: 

Sn _i/2 

— = n ' Poisson 

(1) 

= n -5 / 6 Discreteness 

The former comes from random 
placement of particles; the latter 
from random displacement of par- 
ticles within the cells of a uniform 
mesh. It was known that the first 
gave too large an amplitude, so 
that the initial conditions had to be 
sub-random, but the particle-in-box 
spectrum has too small an ampli- 
tude. There was a guess that some- 
thing like the Zeldovich spectrum 
was needed (Sn/n = 10 — 4 n — 2 / 3 ), 
but there were no ideas on where it 



might come from. 

The situation today is much more 
healthy; 16 years have passed since 
it was realized that inflation could 
seed these initial inhomogeneities, 
and a variety of models of inflation- 
ary structure formation have been 
analyzed in detail. Furthermore, an 
alternative mechanism for seeding 
perturbations exists through the ac- 
tion of topological defects. 

In its simplest form (see e.g. Liddle 
& Lyth 1993; Lyth & Riotto 1998), 
inflation deals with a single quantum 
scalar field, whose expectation value 
rolls down a potential (figure 1). The 
equation of motion is 
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(2) 



If the potential is flat in the sense 
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then vacuum-dominated 'slow-roll' 
inflationary expansion can happen, 
leading to the following horizon-scale 
amplitude for the density fluctua- 
tions: 
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Figure 1. Inflation claims that fluctuations are seeded by quantum fluctua- 
tions in the rolling of a scalar field, <f>, down a potential, V(<j>). 



There are a few free parameters 
here, but the characteristic predic- 
tion of inflation is that there is a 
background of gravity waves. There 
are thus both scalar and tensor con- 
tributions to the CMB power spec- 
trum: 

Cf <x r~ 3 , Cj <x (6) 

(beware: different definitions of the 
tensor index exist). The slopes of 
these spectra depend on the infla- 
tionary parameters: 



scalar tilt = 1 



n$ = 



rflng| 
d In k 



= 6e-2n 
tensor tilt = 1 — n T = 2e 



(7) 



and so does the relative tensor contri- 
bution, giving the consistency equa- 
tion (Starobinsky 1985; Lidsey et al. 
1997): 



ci 



12e = 6(1 - n T ) 



(8) 



These predictions can be violated 
in models with several scalar fields, 
so inflation is a slightly soft target. 
Nevertheless, it would be striking if 
something like these relation was to 
be found in practice. To have both 
a potential mechanism for generat- 
ing the inhomogeneous universe plus 
anything approaching a test of the 
theory is an astonishing achievement 
— whether or not the theory has any- 
thing to do with reality. 

In testing a theory, it is good to 
have an alternative in mind, but I 
will say little about the main class 
of alternatives, which are topologi- 
cal defects. The state of play here 
will be reviewed in Pen's talk, but 



presently looks unpromising (e.g. Al- 
brecht, Battye & Robinson 1998). 



2 The fluctuation spectrum 

2.1 The linear spectrum 

The basic picture of inflationary 
models is therefore of a primordial 
power-law spectrum, written dimcn- 
sionlessly as the logarithmic contri- 
bution to the fractional density vari- 
ance, <t 2 : 



A 2 (fc) 



da 2 
dink 



k i+r 



(9) 



where n stands for ns hereafter. This 
undergoes linear growth 



5 k (a) = <5fc(ao) 



D(a) 



-D(a ) 

where the linear growth law is 
D(a) = ag[Cl(a)] 



(10) 



(11) 



in the matter era, and the growth 
suppression for low S7 is 



g(n) ~ f! 0,65 (open) 
~ n"- 23 (flat) 



(12) 



The transfer function depends on 
the dark-matter content as shown in 
figure 2. 

Note the baryonic oscillations in 
figure 2; these can be significant even 
in CDM-dominatcd models when 
working with high-precision data. 
Eiscnstcin & Hu (1998) arc to be 
congratulated for their impressive 
persistence in finding an accurate 
fitting formula that describes these 
wiggles. This is invaluable for carry- 
ing out a search of a large parameter 



Models for large-scale structure 



3 




Figure 2. Transfer functions for various dark-matter models. The scaling 
with Qh 2 is exact only for the zero-baryon models; the baryon results are 
scaled from the particular case Q B = 1, h = 1/2. 



space. 

The state of the linear-theory 
spectrum after these modifications is 
illustrated in figure 3. The primordial 
power-law spectrum is reduced at 
large k, by an amount that depends 
on both the quantity of dark matter 
and its nature. Generally the bend 
in the spectrum occurs near 1/fc of 
order the horizon size at matter- 
radiation equality, oc (Qh 2 )" 1 . For 
a pure CDM universe, with scale- 
invariant initial fluctuations (n = 1), 
the observed spectrum depends only 
on two parameters. One is the shape 
r = Qh, and the other is a normal- 
ization. On the shape front, a gov- 
ernment health warning is needed, as 
follows. It has been quite common to 
take T-based fits to observations as 
indicating a measurement of Qh, but 
there are three reasons why this may 
give incorrect answers: 

(1) The dark matter may not be 
CDM. An admixture of HDM will 
damp the spectrum more, mimicking 
a lower CDM density. 

(2) Even in a CDM-dominated 
universe, baryons can have a signifi- 
cant effect, making Y lower than Qh. 
An approximate formula for this is 
given in figure 3 (Peacock & Dodds 
1994; Sugiyama 1995). 



(3) The strongest (and most- 
ignored) effect is tilt: if n ^ 1, then 
even in a pure CDM universe a T- 
model fit to the spectrum will give a 
badly incorrect estimate of the den- 
sity (the change in Qh is roughly 
0.3(n - 1); Peacock & Dodds 1994). 

2.2 Normalization 

The other parameter is the nor- 
malization. This can be set at a 
number of points. The COBE nor- 
malization comes from large-angle 
CMB anisotropies, and is sensitive 
to the power spectrum at k ~ 
10~ 3 /iMpc -1 . The alternative is to 
set the normalization near the quasi- 
linear scale, using the abundance 
of rich clusters. Many authors have 
tried this calculation, and there is 
good agreement on the answer: 

as ~ (0.5- 0.6) n~ ' 6 . (13) 

(White, Efstathiou & Frcnk 1993; 
Eke ct al. 1996; Viana & Liddle 
1996). In many ways, this is the most 
sensible normalization to use for LSS 
studies, since it does not rely on an 
extrapolation from larger scales. 

Within the CDM model, it is al- 
ways possible to satisfy both these 
normalization constraints, by appro- 
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COBE ; '"' 




T = Qh exp[-Q B (l + V2h/Q)] 



Apparent density from matter-radiation horizon 



wavenumber k 



a\ = A 2 (A: eff ), k eS / h MpcT 1 = 0.172 + 0.011 [ln(r/0.34)] 2 



Figure 3. This figure illustrates how the primordial power spectrum is mod- 
ified as a function of density in a CDM model. For a given tilt, it is always 
possible to choose a density that satisfies both the COBE and cluster normal- 
izations. 



priate choice of T and n. This is il- 
lustrated in figure 4. Note that vac- 
uum energy affects the answer; for 
reasonable values of h and reasonable 
baryon content, flat models require 
Clm — 0.3, whereas open models re- 
quire S7 m ~ 0.5. 

2.3 The velocity spectrum 

Given the density spectrum, other 
quantities of interest can be calcu- 
lated in linear theory. For example, 
the velocity spectrum is 

A plot of this function is shown 
in figure 5, determined in various 
ways. If these determinations are all 
normalized according to the cluster 
abundance, there appears to be rela- 
tively little uncertainty in the veloc- 
ity spectrum. This is not surprising, 



since the cluster normalization is tied 
to the observed abundance of a set of 
objects of a fixed observed velocity 
dispersion. 

Note the impressive breadth of 
the velocity spectrum: even modes of 
600 ?i — 1 Mpc wavelength contribute 
about 100 kms - 1 rms to the motion 
of the local group. This shows why 
it has been hard to obtain redshift 
surveys of sufficient depth to obtain 
a convergent prediction of the local 
gravitational dipolc. 

These velocity fields arc responsi- 
ble for the distortion of the cluster- 
ing pattern in redshift space, as first 
clearly articulated by Kaiser (1987). 
For a survey that subtends a small 
angle (i.e. in the distant-observer 
approximation), a good approxima- 
tion to the anisotropic rcdshift- 
space Fourier spectrum is given by 
the Kaiser function together with a 
damping term from nonlinear effects: 

5 fe = S r k (l + f3n 2 )D(kan), (15) 
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Figure 4. For 10% baryons, the value of n needed to reconcile COBE and 
the cluster normalization in CDM models. 



where j3 = S7^ 6 /ft, b being the lin- 
ear bias parameter of the galaxies 
under study, \x = k • f , and D(y) ~ 
(1 + y 2 /2)- 1 / 2 , and a ~ 350 kms" 1 
is the pairwise velocity dispersion 
of galaxies (e.g. Ballinger, Peacock 
& Heavens 1996). In principle, this 
distortion should be a robust way 
to determine f2 (or at least 0); see 
Strauss & Willick (1995) and Hamil- 
ton (1997) for details of the practical 
application of these ideas. 

2.4 The nonlinear spectrum 

On smaller scales (fc > 0.1), 
nonlinear effects become important. 
These are relatively well understood 
so far as they affect the power spec- 
trum of the mass (e.g. Hamilton et al. 
1991; Jain, Mo & White 1995; Pea- 
cock & Dodds 1996). These methods 
were applied by Peacock & Dodds 
(1994; PD94) to estimate the lin- 
ear spectrum from the observed spec- 
trum of a number of tracers (figure 
6). To within a scatter of perhaps a 
factor 1.5 in power, the results were 
consistent with a V ~ 0.25 CDM 
model. Even though the subsequent 
sections will discuss some possible 
disagreements with the CDM mod- 
els at a higher level of precision, the 
general existence of CDM-like curva- 
ture in the spectrum is likely to be an 



important clue to the nature of the 
dark matter. 

One interesting aspect of the exist- 
ing spectrum determinations is that 
they are relatively smooth, whereas 
figure 2 shows that rather large oscil- 
latory features would be expected if 
the universe was baryon dominated. 
The lack of such features is one rea- 
son for believing that the universe 
might be dominated by collisionlcss 
nonbaryonic matter (consistent with 
primordial nucleosynthesis if Q m > 
0.1). 

Nevertheless, baryonic fluctua- 
tions in the spectrum can become 
significant for high-precision mea- 
surements. Figure 7 shows that or- 
der 10% modulation of the power 
may be expected in realistic bary- 
onic models (Eisenstein & Hu 1998; 
Goldberg & Strauss 1998). Most of 
these features are however removed 
by nonlinear evolution. The highest- 
k feature to survive is usually the 
second peak, which almost always 
lies near k = 0.05 Mpc -1 (no h, 
for a change). This feature is rel- 
atively narrow, and can serve as a 
clear proof of the past existence of 
baryonic oscillations in forming the 
mass distribution (Mciksin, White & 
Peacock 1998). However, figure 7 em- 
phasises that the easiest way of de- 
tecting the presence of baryons is 



6 



J. A. Peacock 




Figure 5. The velocity power spectrum, according to linear theory. The two 
solid curves show CDM models with Q = 1 and 0.3 (flat), with h = 0.65 
and the appropriate value of n chosen as in figure 4 in order to satisfy both 
the cluster- abundance erg and the CO BE normalization. The points show 
the APM power spectrum linearized as described in Peacock (1997), for the 
same two cosmologies, and also with the cluster normalization. All these 
determinations are in general agreement over the amplitude and broad range 
of the velocity spectrum. 



likely to be through the CMB spec- 
trum. The oscillations have a much 
larger 'visibility' there, because the 
small-scale CMB anisotropics come 
directly from the coupled radiation- 
baryon fluid, rather than the small- 
scale dark matter perturbations. 

2.5 The APM spectrum 

In the past few years, much atten- 
tion has been attracted by the esti- 
mate of the galaxy power spectrum 
from the APM survey (Baugh & Ef- 
stathiou 1993, 1994; Maddox ct al. 
1997). The APM result was gener- 
ated from a catalogue of ~ 10 6 galax- 
ies; because it is based on a deprojec- 
tion of angular clustering, it is im- 
mune to the complicating effects of 
redshift-space distortions. 

The error estimates are derived 
empirically from the scatter between 
independent regions of the sky, and 
so should be realistic. If there are no 
undetected systematics, these error 
bars say that the power is very accu- 



rately determined, and this has sig- 
nificant implications if true. 

A number of authors have pointed 
out that the detailed spectral shape 
appears to be inconsistent with that 
of nonlinear evolution from CDM 
initial conditions, (e.g. Efstathiou, 
Sutherland & Maddox 1990; Klypin, 
Primack & Holtzman 1996; Peacock 
1997). Perhaps the most detailed 
work was carried out by the VIRGO 
consortium, who carried out N = 
256 simulations of a number of 
CDM models (Jenkins et al. 1998). 
Their results are shown in figure 
8, which gives the nonlinear power 
spectrum at various times (cluster 
normalization is chosen for 2 = 0) 
and contrasts this with the APM 
data. The lower small panels are the 
scale-dependent bias that would re- 
quired if the model did in fact de- 
scribe the real universe, defined as 

/A 2 , (fe)\ 1/2 
K*)^ -p • (16) 
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Figure 6. The PD94 compilation of power-spectrum measurements. The up- 
per panel shows raw power measurements; the lower shows these data cor- 
rected for relative bias, nonlinear effects, and redshift- space effects. 



In all cases, the required bias is non- 
monotonic; it rises at k > 5 h^ 1 Mpc, 
but also displays a bump around 
k ~ 0.1 h^ 1 Mpc. If real, this feature 
seems impossible to understand as a 
genuine feature of the mass power 
spectrum; certainly, it is not at a 
scale where the effects of even a large 
baryon fraction would be expected to 
act (Eisenstein et al. 1998; Meiksin, 
White & Peacock 1998). 

An alternative way of presenting 
this difficulty for CDM models is 
shown in figure 9. Here, the meth- 
ods of Peacock & Dodds (1996) have 
been used in an attempt to recover 
the empirical linear spectrum from 
the APM data. It is assumed that 
the bias is independent of scale, and 
that the cluster normalization ap- 



plies. The results (for low and high 
Q) arc contrasted with CDM mod- 
els with reasonable Hubble parame- 
ters and baryon fractions. Interest- 
ingly, there is not a huge difference 
in shape between the CDM mod- 
els. This arises because both are also 
consistent with COBE, and there- 
fore have very different values of the 
tilt, as above. This plot should drive 
home the lesson that the apparent 
value of r cannot be used as an esti- 
mate of Qh without further, possibly 
fragile, assumptions. 

In any case, the linear spectrum 
that would evolve into the APM re- 
sult is seen in figure 9 to differ signif- 
icantly from the CDM models. The 
bump at k ~ 0.1 h^ 1 Mpc is appar- 
ent, as is a flattening at larger k — 
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Figure 7. Baryonic fluctuations in the spectrum can become significant for 
high- precision measurements. Although such features are much less impor- 
tant in the density spectrum than in the CMB (first panel), the order 10% 
modulation of the power is potentially detectable. However, nonlinear evolu- 
tion has the effect of damping all beyond the second peak. This second feature 
is relatively narrow, and can serve as a clear proof of the past existence of 
oscillations in the baryon-photon fluid (Meiksin, White & Peacock 1998). 



the effective spectral index appears 
to approach n = —3, whereas the 
CDM models are never this flat. 

A spectrum of this form would be 
expected with a roughly 30% admix- 
ture of massive neutrinos in a total 
fi ~ 1 (Klypin et al. f993; Pogosyan 
& Starobinsky 1995). This success 
in matching the empirical spectrum 
shape is a very significant attraction 
of the MDM model. However, there 
are also very serious difficulties with 
the idea of an MDM universe, related 
to the formation of high-redshift ob- 
jects. MDM suffers from the same 
difficulties as any Q = 1 model in 
making high-redshift massive clus- 
ters (e.g. Henry 1997). This alone 
may be a strong enough argument 
to reject the model. However, the 
very flatness of the high-fc spectrum 
that is attractive in terms of clus- 
tering statistics makes it difficult for 
MDM to assemble the observed high- 



redshift galaxies (e.g. Mo & Miralda- 
Escude 1994; Mo & Fukugita 1996; 
Peacock et al. 1998). These problems 
are eased in a low-density model, 
which is probably required by the su- 
pernova Hubble diagram (Perlmutter 
et al 1997; Riess et al. 1998), but the 
MDM spectral shape is then wrong. 

3 Towards an understanding 
of bias 

The conclusions from the above 
discussion are either that the physics 
of dark matter and structure forma- 
tion are more complex than in CDM 
models, or that the relation between 
galaxies and the overall matter dis- 
tribution is sufficiently complicated 
that the effective bias is not a simple 
slowly- varying monotonic function of 
position. 



Models for large-scale structure 



9 




Figure 9. a comparison of the linearized APM power spectrum with two CDM 
models that satisfy COBE and cluster normalization. 
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3.1 Local bias 

The simplest assumption is that 
all the complicated physical effects 
leading to galaxy formation depend 
in a causal (but nonlinear) way on 
the local mass density, so that we 
write 

Plight = /(pmass). (17) 

Coles (1993) showed that, under 
rather general assumptions, this 
equation would lead to an effective 
bias that was a monotonic function 
of scale. This issue was investigated 
in some detail by Mann, Peacock & 
Heavens (1998), who verified Coles' 
conclusion in practice for simple few- 
parameter forms for /, and found in 
all cases that the effective bias varied 
rather weakly with scale. The APM 
results thus are either inconsistent 
with a CDM universe, or require non- 
local bias. 

A puzzle with regard to this con- 
clusion is provided by the work of 
Jing, Mo & Borner (1998). They 
evaluated the projected real-space 
correlations for the LCRS survey (see 
figure 10). This statistic also fails to 
match the prediction of CDM mod- 
els, but this can be amended by in- 
troducing a simple antibias scheme, 
in which galaxy formation is sup- 
pressed in the most massive haloes. 
This scheme should in practice be 
very similar to the Mann, Peacock 
& Heavens recipe of a simple weight- 
ing of particles as a function of the 
local density; indeed, the main effect 
is a change of amplitude, rather than 
shape of the correlations. The puzzle 
is this: if the APM power spectrum is 
used to predict the projected corre- 
lation function, the result agrees al- 
most exactly with the LCRS. Either 
projected correlations are a rather 
insensitive statistic, or perhaps the 
Baugh & Efstathiou dcconvolution 
procedure used to get P(k) has ex- 
aggerated the significance of features 
in the spectrum. The LCRS results 
are one reason for treating the appar- 
ent conflict between APM and CDM 
with caution. 

3.2 High-order correlations 

A further piece of evidence con- 
cerning bias comes from the higher- 
order correlations of the density field, 



which are most conveniently de- 
scribed via the hierarchical moments: 

Sj = (5 J )/[(P)] J - 1 . (18) 

These are defined so that, for 
gravitationally-driven growth of fluc- 
tuations, the Sj arc dimensionlcss 
numbers that can be calculated in 
perturbation theory. For example, if 
the density field is smoothed with a 
spherical top-hat window, the third 
moment for a Gaussian field should 
be 

S3 = y - (3 + n) (19) 

(Bernardeau 1994). In the presence 
of pure linear bias of the density field, 
5 g ais = b<5 mass , Sj — ► b 2 ~ J Sj, so 
there might seem to be an interesting 
test here both of Gaussianity and of 
whether light traces mass. In prac- 
tice, Gaztanaga & Fricman (1994) 
found that higher-order results from 
(again) the APM catalogue were in 
remarkably good agreement with the 
simplest unbiased Gaussian predic- 
tions. However, the interpretation of 
this result needs care: linear den- 
sity bias is unphysical, in that it 
can yield 5 < — 1. Gravitational 
evolution causes increasing skewness 
to develop because this is the only 
way in which a large second moment 
(<5 2 ) 3> 1 can coexist with the S > — 1 
constraint. Any bias scheme that has 
the effect of increasing (<5g als ) rela- 
tive to (<5mass) i s almost certain to in- 
crease the skewness also. The higher- 
order results are thus quite a signif- 
icant limit on non-Gaussianity, but 
do not strongly constrain bias. 

3.3 Halo correlations 

In reality, bias is unlikely to be 
completely causal, and this has led 
some workers to explore stochastic 
bias models, in which 

Plight = /(Pmass) + €, (20) 

where e is a random field that is 
uncorrclated with the mass density 
(Pen 1998; Dekcl & Lahav 1998). Al- 
though truly stochastic effects are 
possible in galaxy formation, a rela- 
tion of the above form is expected 
when the galaxy and mass densities 
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Figure 10. The projected correlation function from the LCRS fails to match 
CDM models when comparison is made to just the mass distribution. How- 
ever, the agreement is excellent when allowance is made for a small degree of 
scale- dependent antibias; galaxy formation is suppressed in the most massive 
haloes (Jing, Mo & Borner li 



are filtered on some scale (as they al- 
ways are, in practice). Just averaging 
a galaxy density that is a nonlinear 
function of the mass will lead to some 
scatter when comparing with the av- 
eraged mass field; a scatter will also 
arise when the relation between mass 
and light is non-local, however, and 
this may be the dominant effect. 

The simplest and most impor- 
tant example of non-locality in the 
galaxy-formation process is to rec- 
ognize that galaxies will generally 
form where there are galaxy-scale 
haloes of dark matter. In the past, it 
was generally believed that dissipa- 
tive processes were critically involved 
in galaxy formation, since pure col- 
lisionless evolution would lead to 
the destruction of galaxy-scale haloes 
when they are absorbed into the cre- 
ation of a larger-scale nonlinear sys- 
tem such as a group or cluster. How- 
ever, it turns out that this overmerg- 
ing problem was only an artefact of 
inadequate resolution. When a sim- 
ulation is carried out with ~ 10 6 
particles in a rich cluster, the cores 



of galaxy-scale haloes can still be 
identified after many crossing times 
(Ghigna et al. 1997). This is a very 
important result, since it holds out 
the hope that many of the issues con- 
cerning where galaxies form in the 
cosmic density field can be settled 
within the domain of collisionless 
simulations. Dissipative physics will 
still be needed to understand in de- 
tail the star-formation history within 
a galaxy-scale halo. Nevertheless, the 
idea that there may be a one-to- 
one correspondence between galaxies 
and galaxy-scale dark-matter haloes 
is clearly an enormous simplification 
— and one that increases the chance 
of making robust predictions of the 
statistical properties of the galaxy 
population. 

Another piece of work which also 
emphasises the survival of galaxy- 
scale haloes is Klypin et al. (1997; 
see also Kravtsov's talk at this meet- 
ing). Rather than looking only at a 
single cluster, Klypin et al. attempt 
to find galaxy haloes within a ran- 
dom cosmological volume. Their re- 
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Figure 11. Klypin et al. (1997) have shown that the correlation function of 
galaxy-scale haloes is much closer to a single power law than the correlation 
function of the overall mass distribution. 



suits arc shown in figure 11, where 
the correlation function for the mass 
in a ACDM simulation is contrasted 
with the correlation function for the 
haloes. The mass £(r) shows the fa- 
miliar form: the sharp quasilinear 
rise between £ ~ 10 and 300, followed 
by the flattening associated with viri- 
alization. 

By contrast, the halo correlation 
functions are much closer to single 
power laws: they are antibiased at 
£ ~ 100, but positively biased for 
£ > 1000. These are qualitatively the 
trends that are required in order to 
reconcile the APM data with CDM 
models, so this is clearly a result of 
great potential significance. Future 
work should improve the precision 
of the results and allow the simula- 
tion boxes to be made larger (they 
are presently uncomfortably small). 
However, if the result holds up, we 
will be no nearer an understanding 
of what is going on. A simple re- 
sult like power-law correlations de- 
mands to be explained in a more in- 



tuitive manner. If the result turns 
out to be inevitable on some general 
grounds, then this may be bad news 
for the idea of large-scale-structure 
studies as a diagnostic tool - it is 
possible that the galaxy correlations 
may turn out to be less sensitive 
to the underlying linear power spec- 
trum than the mass correlations are. 



4 Clustering at high redshifts 

A longstanding ambition has been 
to remove some of the uncertain- 
ties concerning the generation of 
large-scale structure by observing 
the growth of clustering with time. 
After a number of years of ambigu- 
ous results from angular clustering of 
faint galaxies, this has at last proved 
possible by using deep redshift sur- 
veys. 

Perhaps the first convincing mea- 
surement of clustering evolution 
came from the CFRS (Le Fevre et 
al. 1996), which indicated evolution 
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in the sense that £ at a given co- 
moving separation was a factor of 
about four smaller at z = 1 than it 
is today. Evolution in the same sense 
(but somewhat less strong) has been 
found by the CNOC2 team (Carlberg 
ct al. 1998). 

These results appear to prove that 
clustering does indeed grow by grav- 
ity, but this idea looks less con- 
vincingly demonstrated when the re- 
sults from clustering at higher red- 
shifts are considered. It is remark- 
able that samples at z ~ 3.5 are now 
as large as those at 0.5 % z < 1 
were only two years ago; this shows 
how efficient the Lyman-break se- 
lection technique has been. At this 
highest rcdshift, we have the remark- 
able result that the comoving corre- 
lation strength appears to be very 
close to that measured today (Steidel 
ct al. 1998; Adelberger et al. 1998). 
In other words, any downward evo- 
lution in £ between z = and 2 = 1 
must be reversed between z = 1 and 
z = 3.5, so that £(2) has a minimum. 

Alternatively, if we believe that 
clustering of the mass grows mono- 
tonically according to gravitational 
instability, then clustering at 2 = 3.5 
must be very strongly biased. Stei- 
del et al. (1998) obtain a minimum 
bias of 6, with a preferred value of 8, 
assuming Q = 1. It is easy to trans- 
late to other models, since we observe 
cell variances c 2 ^ directly, where the 
cell has a given angular width and 
depth in rcdshift. For low Q mod- 
els, the cell volume will increase by 
a factor [5|(r) dr\/[S 2 (ri) dn]; com- 
paring with present-day fluctuations 
on this larger scale will tend to in- 
crease the bias. However, for low fi, 
two other effects increase the pre- 
dicted density fluctuation at z = 3: 
the cluster constraint increases the 
present-day fluctuation by a factor 
SI -0,56 , and the growth between red- 
shift 3 and the present will be less 
than a factor of 4. Applying these 
corrections gives 

*»(" = 0-3) = f 0.42 (open) 
b(U = 1) 1 0.60 (flat) ' 1 ' 

which suggests an approximate scal- 
ing as b oc n - 72 (open) or U 0A2 
(fiat). Strong bias is needed at 2 = 3 
for all reasonable values of Q. 

The standard explanation for this 
effect is high-peak bias, which is 



bound to operate if the high-redshift 
galaxies are rare, newly-forming sys- 
tems. The linear bias parameter de- 
pends on the rareness of the fluctu- 
ation and the rms of the underlying 
field as 

v 1 - 1 v 2 - 1 
6 = 1+ = 1 + — (22) 

1/(7 C 

(Kaiser 1984; Cole & Kaiser 1989; 
Mo & White 1996 - but sec also 
Jing 1998), where v = 8 c /cr, and 
a 2 is the fractional mass variance 
at the redshift of interest. A bias 
at the observed level is therefore to 
be expected, provided the Lyman- 
break galaxies are moderately mas- 
sive systems, with circular velocities 
V ~ 150 kms- 1 (e.g. Bagla 1998a, b; 
Peacock et al. 1998). The minimum 
in § makes sense since we observe 
something proportional to (ba) 2 ; b 
is constant for large <r, but scales as 
<7~ 2 for small <r; the observed cluster- 
ing therefore has a minimum at the 
rcdshift where <r ~ 1. 

In order to turn this plausibil- 
ity argument into a proof, two tests 
should be carried out. First, the 
circular velocities of Lyman-break 
galaxies need to be measured; sec- 
ond, even larger and more detailed 
rcdshift surveys need to be carried 
out at 2 ~ 3, in order to verify that 
the structures causing £(r) have the 
same topological character of voids 
and filaments as seen in local large- 
scale structure. If these criteria are 
satisfied, then we can claim to have 
some understanding of the growth of 
structure, and the existing clustering 
results may be used as a test of mod- 
els. 



5 Conclusions 

It should be clear from this talk 
that large-scale structure has ad- 
vanced enormously as a field in the 
past two decades. Many of our long- 
standing ambitions have been re- 
alised; in some cases, much faster 
than we might have expected. Of 
course, solutions for old problems 
generate new difficulties. Although 
we now have good measurements of 
the clustering spectrum and its evo- 
lution, it is less clear that the place- 
ment of galaxies with respect to the 
mass is simple enough to use these 
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results as direct statistics with which 
to test theories. A fairly safe bet is 
that one of the major results from 
new large surveys such as 2dF and 
Sloan will be a heightened apprecia- 
tion of the subtleties of this problem. 

Nevertheless, we should not be de- 
pressed that problems remain. Ob- 
servationally, we are moving from an 
era of 20% - 50% accuracy in mea- 
sures of large-scale structure to a fu- 
ture of pinpoint precision. This ma- 
turing of the subject will demand 
more careful analysis and rejection of 
some of our existing tools and habits 
of working. The prize for rising to 
this challenge will be the ability to 
claim a real understanding of the 
origin of structure in the universe. 
We are not there yet, but there is 
a real prospect that the next 5—10 
years may see this remarkable goal 
achieved. 
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